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A general Reynolds analogy (GRA) theory is proposed for the mean and fluctuating 
velocity and temperature in compressible wall-bounded turbulent flows. In particular, 
an exact analogy solution is derived for compressible turbulent pipe and channel flows 
and an approximate analogy solution is derived for compressible turbulent boundary 
layers (CTBL), both of which are independent of fluid Prandtl number and wall temper- 
ature condition. The analogy solutions are in excellent agreement with direct numerical 
simulation data, able to reproduce empirical relations, and can be viewed as exten- 
sions of existing theories. In contrast to Walz's equation for adiabatic CTBL, the mean 
temperature-velocity relation derived by GRA can be applied to different wall-bounded 
flows in non- adiabatic wall condition, which is achieved by extending Walz's adiabatic 
recovery factor to a heat flux dependent one. The fluctuation temperature-velocity rela- 
tions derived by GRA are slightly different from the modified strong Reynolds analogy 
derived phenomenologically by Huang et al. (HSRA), and have a better performance than 
HSRA. In addition, several key quantities are introduced in GRA, including a general 
total enthalpy (or temperature) and an adiabatic degree-a well-founded dimensionless 
parameter for characterizing the wall-temperature effects in non-adiabatic flows. The 
GRA unveils the universal feature behind the complex nonlinear couplings between the 
thermal and velocity fields, and makes possible of predicting the mean fields of compress- 
ible wall-bounded turbulence with the information of the corresponding incompressible 
flow. 

Key words: Compressible turbulence. Shear layer turbulence. Turbulence theory, Com- 
pressible boundary layers 



1. Introduction 

There are strong nonlinear couplings between the velocity field and the temperature 
field of compressible wall-bounded turbulent flows (jGaviglid Il987t ISmits fc Dussauge 
The temperature-velocity relationship is crucial to the design of hig h speed ve- 



hicles; thus it has been the focus of many studies during the past decades ( "Wald 19691 : 
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theless, it remains an unresolved question owing to the complex character of compress- 
ible turbulence, and owing to the various wall temperature conditions (adiabatic/non- 
adiabatic walls) and flow situations (external/internal flows) involved in aerospace engi- 
neering. 

The first theory on the temperature- velocity relation ship was present ed by O. Reynolds 
in 1874 for incompressible laminar boundary layers ( Revnolda 1961 ). Then it was ex- 
tended to incompressible turbulent boundary layers (ITBL). For a zero-pressure-gradient 
(ZPG) thermal ITBL, the extension gives a linear dependence of temperature on velocity, 
a result of the similarity betwee n the momentu m equation and the energy equation, i.e. 
the so-called Reynolds analogy (|Gavigliol 119871) . The similarity is largely attr ibuted to 
the negligible viscous dissipation in incompressible turbulence ( Gavigliol 1987 ). In com- 
pressible flow, however, the viscous dissipation plays an important role (|Gaviglid 11987 : 
ISmits fc Dussauge.,2006i) . prohibiting a straightforward analogy. 

A lot of studies were conducted to extend the Reynolds analogy to compressible flows. 
In terms of the relationship between mean temperature and mean velocity, several rela- 
tions were derived based o n various approximat ions. Busemann (1931) and Crocco (1932) 
obtained the first relation ([Pirozzoli et a/.ll2004l ) for ZPG compressible laminar boundary 
layers by assuming unity Prandtl number citemorkovinl962effects,gatski2009compressibility 
(for air, Pr w 0.7). Their derivation can be applied to turbulent flows. In contrast to the 
linear relation of Reynolds analogy, the Crocco-Busemann relation exhibits a quadratic 
dependence of mean temperature on mean velocity: 
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where u is streamwise velocity, T is temperature, Cp is specific heat at constant pressure, 
— denotes Reynolds average and f is called extended recovery factor here to distinguish 
from the classical recovery factor r( (Wa lz 1969). Throughout this paper, subscript w 
denotes wall and 5 denotes boundary layer thickness or channel/pipe center. In the 
Crocco-Busemann relation, f is unity, which was modified to r « 0.9 by Walz (1966) to 
account for the deviation of from unity (Walz 1969). Duan et al. (jPuan fc Martinll201ll) 
recently presented an empirical relation based on a vast amount of direct numerical 
simulation (DNS) data up to hypersonic regime. Duan's relation is in excellent agreement 
with DNS data and is valid for various flows, such as flows with low or high enthalpy 
and with or without chemical reaction. The only difference between Walz's equation and 
Duan's relation is that the first u/us in Eq. (|l.ip is replaced by a quadratic function 
of u/us- As to be demonstrated in this paper, Duan's relation can be derived from our 
new theory, called general Reynolds analogy (GRA). The GRA gives a general mean 
temperature-velocity relation with a variable f that is only associated with the wall- 
temperature condition and can be expr essed as a functi on of r, Pr, s and 0, where s is 
the well-known Reynolds analogy factor ( Bradshawll977t ) and 9 — {Tw—Ts)/{Tr—Ts) is a 
newly defined adiabatic degree based on the recovery temperature (i.e. the adiabatic wall 
temperature) Tr — Ts + rug/{2Cp). For adiabatic wal ls, the gene ral mean temperature- 
velocity relation naturally reduces to Walz's equation! Walj 19691 ) . 

As to the fiuctuation fields, a series of relations between the streamwise velocity fluc- 
tuation u and the temperature fluctuation T were also derived. The first relations were 
identified by Morkovin in 1962 (primarily due to Young in 1951 ( Spina et al\\l99^ ) and 
were known collectively as Strong Reynolds Analogy ( SRA). SRA w ere derived for ZPG, 
adiabatic CTBL under two assumptions: (a) Pr — 1 (jGuarini et a/..,200Q) or Prm = 1 
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. . , ______ ,_, _^ (|Morkovinlll962l ) or 

200Cl[ ). where Tf = T + UiUi/(2Cp) is total temperature. In 1974 SRA was extended 
to non-adiabatic walls by Cebeci and Smith (|Cebeci fc SmithI 1974) and was c alled ex- 
tended SRA (ESRA) Owing to the too rough assumptions (Morkovin"l962'; 'Gavigliql 
1987tlSpina a?l ll994 iGuarini et aZll2000l : iMaeder et alll2001; Smits fc Dussauge 2006i: 



(Guarini et al 



Duan et al. 2010l ). bo th SRA and ESRA deviate clearly from experimental (iGavigh 



198d) and DNS data (jGuarini et aLll2000l: iMaeder et all 120011 : iDuan aLll201C 
hence are inappropriate for applying to turbulence models. To avoid these inappropriate 
assumptions, some modified SRA were proposed by several authors, including Gaviglio 
(GSRA,1987), Rubesin(RSRA,1990) and Huang (HSRA,1995). These modified SRA have 
a common form; 



(7- 1) A/Vu'Vu 
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where c = 1 for GSRA, c = 1.34 for RSRA and c = Pr* for HSRA, 7 = Cp/Cy = 1.4 
is the ratio of specific heat, M — uj \J jRT is local Mach number. Among the three 
relations, HSRA best agrees with DNS data at various wall temperature conditions and 
for various flows includ ing C TBL ICuarini a^.ll2q00t IMaeder et all l200lt iPirozzoli et al\ 



2004lDuan et aLll2010ll201ll:lDuan fc Marti'nll201l[) and com pressible channel flow (CCF) 
( Huang et al. 19951 ). Guarini et al. ( Guarini et al. 2000l) later presented a thorough 
analysis and pointed out that HSRA revealed the key analogy between the normal- 
ized rates of turbulent heat and momentum transfer. More specifically, the correlation 
coefficients Rv'u' and Rv'T' are nearly equal throughout the boundary layer, where 



Ra' p' = \l ot P' / (v a'^ y 13''^). In this paper, we will prove that this key analogy is a 
consequence of GRA, and a modified HSRA (MHSRA) is resulted, which demonstrates 
a better performance than HSRA. 

All the previous analogy theories hold only for ZPG flows. For flow driven by pres- 
sure gradient, such as channel and pipe flows, the streamwise pressure gradient excludes 
a straightforward application of the analogy theory. As to be shown, GRA presents 
as a unified theory for a resolution of the problem. In addition, even for ZPG CTBL, 
the existing theories have more or less defects. For example, both SRA and the mod- 
ified SRA yield that Ru'T' = —1 for CTBL and Ru'T' = 1 for CCF and compress- 



ible pipe flow (CPF), w hich, however, disagree with experimental ( Guarini et al\ 



2000; 



2OOOI: IMaeder et al 



Smits fc Dussauge 2006) a nd DNS data (iHuang et a/]|l995 : Guarini et al _ _ _ 

200ll:|Pirozzoh et aL .20041: iDuan et aLll2010ll201ll: iDuan fc Martfnll2oTTI ). The disagree 
ment has not been well understood, but can be explained by GRA. 

In this paper, we present the GRA for compressible wall-bounded turbulent flows. An 
exact analogy solution for CCF/CPF and an approximate analogy solution for CTBL 
are obtained, which are consistent with the existing theories and the empirical rela- 
tions, and in excellent agreement with DNS data. For the DNS data, a notation like 
CCFM4.50ei.0Q is used to denote a CCF with M = 4.50 and 9 = 1.00. x,y,z (or 
Xi,X2tX3) are used to denote streamwise (axial), wall-normal (radical) and spanwise (az- 
imuthal) directions in Cartesian (cylindrical) coordinate, and the corresponding velocity 
components are u, v, w (or ui, U2, U3). Throughout the paper, Reynolds decomposition is 
used only, but all the der ivations and conclusions can be extended to Favre decomposition 
( Gatski fc Bonnet I2OO9I) m a straightforward way. This paper is structured as follows. 
A brief review and discussion of the existing Reynolds analogy theories are presented in 
section [2] The GRA is presented in section [Sj followed by the mean temperature- velocity 
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relationship in sec|4] and the fluctuation temperature-velocity relationship in seclSl A 
discussion of GRA is presented in secEl And finally, conclusions are drawn in secEl 



2. A review of various Reynolds analogies 

The time-averaged Navier-Stokes equations read ( Gatski fc Bonnet 20091 ) 



9x,(pul)=0 (2.1) 
dxj (puiUj) = -dxiP + (t^~^7) (2.2) 
dxj (phtUj) = {uiTx^x, ) - dxj iWj) (2.3) 

,where Eq. (|2.ip is continue equation, Eq. (|2.2p is momentum equation and Eq. (|2.3|) is 
energy equation, p is density, p is pressure, ht = C'pTt = h + UiUi/2 is total enthalpy 
and h = CpT is enthalpy, TxiXj—'^P-{Sij — Skk^ij/3) is viscous stress tensor with 5*^^ = 
{dxjUi -\- dxiUj)/2, Qxi — —kdxiT is heat flux, /i is viscosity that is related with thermal 
conductivity k by Pr — Cpp/k. Throughout the paper, the following approximations are 
frequently adopted, which are highly accurate and well accepted: 

UiUi u^^uu u^; — « 2uu ; [pv)' [u^ /2) k, u[pv)' u . (2-4) 

In the following, we briefly review the existing Reynolds analogy theories. To be concise, 
discussions are limited to ZPG turbulent boundary layers. 

2.1. Reynolds analogy for incompressible turbulent boundary layers 

For ZPG incompressi ble turbulent b oundary layers, the momentum and energy equations 
can be simplified to (jGavig liol ll987h 



pudxU + pvdyU = dy{p,dyU — pu'v') (2-5) 

Pudxf + pvdyf = dyifl/ Prdyf) + dy{~pT'v') + $/Cp (2.6) 



where $ = 2fiSijSji is viscous dissipation and negligible in incompressible flows (|Cebeci fc Bradshaw 



19841 ). Above the viscous sublayer, turbulent transport of momentum and heat predom- 



inates over molecular transport, so the molecular transport terms can also be neglected. 
One then obtains: 



pu*dxU* + pv*dyU* — —dypu*'v*' (2.7) 
pu*dxf* + pv*dyf* = ~dypT*'v*' (2.8) 

,where u* — Ui/us and T* — {T ^T^)/ [T^—Tw)- The corresponding boundary conditions 
are (a). = : li* = T* = (b). y = (5 : It* = T* = 1. The similarity between Eq.jlH) 
and Ea. (|2.7p means u* —T*. In the dimensional form, the temperature- velocity relations 
are 

{T-T^)l{Ts^T^)^u/us (2.9) 
t' /{Ts~T^)^u /us. (2.10) 

Eq. (|2.9p and Eq. (|2.10p show a linear dependence of temperature on velocity for both 



the m ean and fluctuation fields. As pointed out by Cebeci and Bradshaw IjCebeci fc Bradshaw 



19841 ). the instantaneous fiuctuation relation of Eq. (|2.10p denotes an 'exact analog; 



:y 

that only applies to an ideal, i.e. unrealistic, situation. A weak form can be derived from 
Eq. dSini): 



T'7(T5 - T™) « Vu'Vui-, (2.11) 
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which is less questionable. 

2.2. Crocco-Busemann relation (1931,1932), SRA (1962), and ESRA (1974) 

In compressible flows, the viscous dissipation can not be neglected, so the 'Reynolds 
analo gy' for incompressible flows cannot be applied directly. Owing to Young ( Howarth 



1953 ) , the energy equation can be written in the form of total enthalpy (jGatski fc Bonnet 



2009 k 



pudxh + pvdyht = dy[{fi/Pr)dyht - ph'^v'] 

+dy[p{l - l/Pr)9j,(S,S,/2 -<</2)]. (2.12) 
Assuming Pr = 1, momentum equation and energy equation show a similarity: 



pudxU + pvdyU — dy{p.dyU — pu'v') (2-13) 

pudxht + pvdyht — dy{fldyht — ph[v'). (2-14) 

To find an analogy solution of Eq. (|2.13p and Eq. (|2.14p . a rather rigorous assumption 
is introduced: hf = Uqu , where is a proportional constant with dimension of velocity. 
Subtracting Eq.^J^ from Eq.^J^ gives 

{pudx + pvdy - dypdy)(ht - Uqu) = 0. (2-15) 

A possible analogy solution of Eq. (j2.15p is ht — Uqu = const in the whole x — y plane. 
Here, the emphasis of 'possible' is because this solution is a sufficient but unnecessary con- 
dition of Eg. (12. 151) . Applying the boundary conditions, one finds Uq = {hts — htw)/us = 
— Pigyw/Tw Then the mean temperature- velocity relation can be derived as: 

f = | + ^^^ + ^(^)^T. = T. + |-. (2.16) 

Ts Ts Ts us Ts \us J 2Cp 

Eq. (|2.16p is called Crocco-Busemann relation since it is similar to the relation derived 
by Crocco (1932) and Busemann (1931) for compressible laminar boundary layers. 

For an adiabatic wall, qyw = gives Uq — 0. The analogy solutions become Tt = 
T,„ and = 0. The former is moderately confirmed by DNS data ljGuarini et al. 



Pirozzoli et al\\2004 ). which show the deviation of the mean total temperature from a 



constant is less than about 7%. The latter leads to: 

Tt ^CpT +uu' kO. (2.17) 
, where the approximations in Eq. (|2.4I) are used. This instantaneous relation has sev- 



eral statistical consequences: u'T' « —{u/Cp)u'^, V T'^ w {u/Cp)Vu"^, and pT'v' 
— {u/Cp)pu'v', and further yields: 



(7 - 1) A/Vm'Vu 



(2.18) 



Ru'T' = « -1 (2.19) 



/rp'2 



Prt = ^"';'(^^/^^) ^ 1 (2.20) 
pT'v' (du/dy) 
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The above relations are derived by using the approximation Tt ~ CpT + v? /2 and 
it's the normal invariance of Tt, saying ft = = fg. Eq. p.l7p - Eq. p.21|) were 
first identified by Morkovin in 1962 and collectively known as stro ng Reynolds analogy . 
Here the use of 'strong' may come from the assumption = Uqu (iGuarini et~al. 2000) 
(or negligible total temperature fluctuation (|Morkovin 1962 : Lele 19941 )). which is too 
rigorous to meet in a real turbulence. ^^^^^^^^^ 

The assumption was later invalidated by DNS data (jGuarini et "aLl l2000l:[Maeder et al 



200 it iDuan et al. 2011), which show that the magnitude of yT^'^ is comparable to 



V T'2 (see Eg. (12. 171) ). Consequently, the SRA relations are poorly satisfied. For example, 
the predictions R^'x' — ~^ and Prt — 1 clearly deviate from the observed values of 
—R,,'^' = 0.5 ^ 0.7 and Prt = 0.7 ^ . 9 in adiabatic CTBL up to hypersonic regime 
(|Guarini et aLll200d iMaeder alll200ll : IPirozzofi 6^011 12004': 'Puan 6^ 0^.1120111 ). How- 
ever, Eq. (12 171) is well satisfied for both experiments (Gaviglio 198'^ and numerical 
simulations ( Guarini et al. I I2OOOI) . Then a question arises as to how Eq. (|2.17p can be sat- 
isfied under an incorrect assumption. As to this, both Debieve (1976) and Gaviglio (1987) 
(|Gavigliolll987l ) demonstrated that Eq. (|2Tf|) was a sufficient but unnecessary condition 
of Eq. ()2.18|) . By rearranging the definition of the total temperature fluctuations as 



* = 1 (2.22) 

(7- l)Af2Vu'VM 

, Guarini (2000) (jGuarini et all I2OOOI ) pointed out that the validity of Eq. (|2T7)) came 
from -^T^ ~ 2T[T, instead of 7^ « 0. Gaviglio (|Gavigliolll987l ) further observed 
that the total temperature fluctuation could be arranged in a more general form as 



1/2 



(2.23) 



Applying Eq STM to Eg. ST^. one obtains i?„,T' = T'^^ /2T"^ - l (|Gavigliol [l987l : 
Guarini et al. Ipirozzoli et a/.ll2004 [). which has a better performance than Eq. ()2.19|) . 

SRA can be extended to non-adiabatic wall flows. Evaluating ht — Uqu = hw at 5, one 
finds Tts — T-uj = us{Uq/Cp). Substituting the relation into the assumption CpT^ ~ Uqu , 
with the aid of the approximations in Eq. (|2.4|) . one obtains the following relations: 



T /T 



(7 - 1) M^u' /u 



1-a 



i-a 



Tts — 

Tts — 
ugu 



(2.24) 



(2.25) 



(7- l)APVu'^/u 

The i nstantaneous relationship of Eq. (|2.24p was proposed by Cebeci and Smith (ICebeci &: Smith 
19741) in 1974. The r.m.s. form of Eq. (P?^ . i.e. Eg. (1^7251) . was suggested and referred 
to as extended SRA (ESRA) by Gaviglio (iGavigliolllbsTf I in 1987. In case of having heat 
fiux at the wall, ESRA has a noteworthy improvement comparing with SRA, but has 
also an observable deviation from real turbulcncc(Gavigho 1987.). 

2.3. Walz's eguation (1966) 

The deviation of Pr from unity is the major reason responsible for the difference be- 
tween Crocco-Busemann relation and DNS data, as shown in fig. ([2]). To improve, Walz 
presented an approximate solution of the Reynolds averaged Navier-Stokes (RANS) equa- 
tions with the assumption of constant mixing Prandtl number, i.e. Prm — const, which 
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is a reasonable approximation for CTBL and other shear flows (jSmits fc Dussaug 
Prm is defined by Pr™ = Cp(/^ 
viscosity and kt = 



Ht)/{k + kt), where fit = {— pu' v') / {du / dy) is eddy 



-CppT") / [dT / dy) is eddy thermal conductivity. 
In Walz's derivation, temperature is assumed a function of u, and RANS equations are 
written in terms of the independent variables x and u, inst ead of a: and ?/. After neglecting: 
most streamwise derivatives, the energy equation reads ( Walz 19691 : Smits fc Dussauge 
200d) 



f[du{duf/ Prrn) + 1/Cp] + (1/Pr„, - l)aaTaQf = 



(2.26) 



, where r is total shear stress. The boundary conditions are [a), u — Q : T — T^, {df / dy)y=o 
dp/dx — (b). u = Us : T = Ts,f = fs = 0- The assumption Prm = c onst makes 
Eq. (|2.26p can be integrated two times (IWalj 119691: ISmits fc Dussaugdl2006l) . Then one 
is able to obtain: 



Tr 



Ts 



-A 



Ts — Tr 
fs 



f,,fr = f6 + r^ 



(2.27) 



where /i, /2 and r are functions of t/tu,, Pr^n and Wal3ll969l) . Using a linear ap- 
proximatio n f/fw ~ 1 — y/S and Pr^ « 0.86 averaged over the boundary layer thickness, 
Walz(Wal^ ll969l) found that /i w u/us, A ~ {u/usf and r w 0.9 ( 0.88 in reference 
(fWalj 119691) '). Then Eq. ^^Tm has the same form as that of Eq. (|l.ip . and was called 
modified Crocco-Busemann relation or Walz's equation. 



2.4. GSRA(1987),RSRA(1990),HSRA(1995) 

To avoid the inappropriate assumptions made in SRA, Gaviglio attempted to establish 
a relationship between u and T based on the recognition of the universal a nd domi- 
nant role of the large scale motion in compressible wall-bounded turbulence. ( Gavigliol 
19871) . He found that the fluctuating velocity and temperature induced by the large 



scale movements are proportional to the local gradient of mean velocity and temper- 
ature, and the corresponding proportional constant are assumed as velocity mixing 

length (.u = \/u^ / {du/ dy) and temperature mixing length It = Vt^ / (du/ dy) , say- 
ing u — £u{du/dy) and T = iridf /dy). In instantaneous and r.m.s. forms, one obtains 

cT /dyf = u /dyU and cV T'^/dyf — V u'^ /dyU, where c = iu/^T- The latter relation, to- 
gether with the approxi mations in Eq. (j2.4|) . yields the modified SRA, i.e. Eq. (|1.2|) . Based 
on experimental data ( Gavigliolll987t ). Gaviglio further assumed c = 1, which gave the 
GSRA.In 1990, Rubesin independently gave an equivalent forni of GSRA with c = 1.34, a 
constant calculated by Huang et.al in 1995 ( Huang et al. 19951) using various turbulence 
modeling constants given by Rubsein (see E g (20) in reference (lRubesinlll990[ )'). In Con- 
trast to the constant value of c, Huang et.al ( Huang et aZ.lll995l ) (1995) showed that c is 



equal to Prt that varies along the wall-normal direction. Indeed, multiplying both sides of 
cT /dyf = u /dyU by pv and averaging, one obtains c = {pv'u'dyf)/{pv'T'dyu) — Prt, 
which gives HSRA. 

Possibly owing to the universality of the large scale movements in wall-bounded turbu- 
lence, HSRA works well for various wall-bounded tu rbulent flows at var ious wall condi- 
tions. In channel flow with cold wall, Huang et.al ( Huang et al. 1995 h 
had an evident improvement over ESRA. In ZPG CTBL, Duan's 



showed HSRA 
DNS data up to 



M 



12 showed that HSRA was effective for adiabatic and non-adiabatic, ca t alytic 



and non- reacting walls at low- and high- enthalpy conditions ( Duan et al. 20 id . 2011 



Duan fc Martfn..2011,) . 
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2.5. The empirical temperature-velocity relation (2011) 

Recently, Duan et.al ( Duan fc Martinlboill ) successfully removed the explicit dependence 
of the temperature-velocity relationship on thermal and chemical models by an empiri- 
cal relation, in which a dimensionless 'recovery enthalpy'is introduced. For a calorically 
perfect gas, ht = CpT, the recovery enthalpy degenerates to a dimensionless recovery 
temperature f;^ defined as T*y = {fry - f^)/{frs - f^) with Try ^ f + ru^/{2Cp). 
Plotting 7"*^^ versus u/us, Duan et.al found that all DNS data collapsed together, inde- 
pendent on freestream Mach number, wall temperature, surface catalysis and enthalpy 
conditions. The best fitting of fry = f{u/us) gives 



/(- 



us 



0.8259. 



us 



us 



(2.28) 



By expanding T*y — f{u/us), Duan et.al (jDuan fc Martini 120111 ) obtained a mean 
velocity -temperature relation in the following form: 



fs fs 



fs 



Ts - Tr 
fs 



u , 

— I ,Tr 

US, 



Ts+r^- (2.29) 



Comparing Eq. (|2.29l) with Walz's eqaution Eq. (|2.27p . the only difference is the replace- 
ment of u/us in the second term of Walz's relation with f{u/us). 



3. General Reynolds analogy for compressible wall-bounded 
turbulence 

Before proceeding, we recall the key approximations and assumptions involved in the 
previous analogy theories, which result in the final differences between the theories and 
real turbulence. In Crocco-Busemann relation, SRA, and ESRA, an approximation of 
Pr = 1 and an assumption oi H — Uqu are introduced. Both of them are not good 
descriptions of real turbulence. Young (1986) pointed out that the difference caused 
by Pr— 1 7^ represents a measurement of the dissimilarity of the two different modes 
of transfer for vector pu and sc alar Tt which do not respond in the same manner to 
change in density and pressure iGavialid {l98^ ). The assumption H = Uqu is also 



an 'exact' analogy similar to Reynolds analogy. A consequence of this 'exact' equation 
of fluctuations, as well as the fluctuation assumptions to derive GSRA and HSRA, is 
the rigorous relation of |i?„'T'| — 1, which is unfortunately not satisfied in a real wall- 
bounded turbulence. As for Wa lz's equation, a lot of DNS data show that it is satisfactory 



in quasi-adia batic CT BL (P irozzoli aLM2004t iDuan et al\\20lu) . but unsatisfactory 
(Duan et aLll20l5lDmin fc Martinll201l[ ) in non-adiabatic flows, as shown in Fig. ([2]). Two 
key approximations are made to derive Walz's equation: (a) Pr™ = const, (b) f/fu, « 
1 — y/S. Since /i, /2 and r are complex functions of Prm and f/f^,, it is difficult to 
quantify the error introduced by these two assumptions, especially for (b). However, our 
DNS data show that 1 — f/f^ w for quasi-adiabatic wall, where n « 1.35 at 

M = 2.25 and n w 1.10 at M = 6.00. For an over cooled wall, a near-wall overshoot of 
the total stress ( f/f^ > 1) is observed, also inconsistent with the linear approximation 
of(b).^ 

All in all, an improved theory should avoid, if possible, those inappropriate approxi- 
mations or assumptions. Besides, the theory should, ideally, yield the empirical relations 
of Duan et al. and Huang et al. In this section, we present a general analogy between the 
Reynolds-averaged momentum and energy equations that meets the above requirements. 
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3.1. Compressible pipe flow 
For a compressible pipe flow, the governing equations of a thermally and calo rically 



perfect gas in cylindrical coordinates are referred to reference ( Ghosh et al\\20ld[ ). Con- 



yentionally, the radial coordinate and velocity are denoted by r and u^, respectively, the 
former of which is used already by recovery factor in this paper. To avoid ambiguity, 
we also use y and v to denote them. This is helpful to demonstrate the similarities and 
differences between pipe and channel flows. Note that some y in this subsection have 
been marked as y (i.e. y — y), which is used to simplify the derivation of GRA in channel 
flow (see subsection. 13. 2p . Then the momentum and energy equations averaged in time 
and in axial and azimuthal directions are written as (see Appendix A): 



ydj:p = dy[y{T^y - {pv)'u')] (3.1) 



^ dy[y{uTx:y ~ (pv)' h[ + vTyy + kdyT)]. (3.2) 



Using approximations of UTxy ~ uTxy, Txy ~ fJ,dyU, Qy = —kdyT !v —kdyT and approx- 
imations in Eq. (|2.4I) . neglecting the small term dy{vTyy) and defining total stress as 

f = jldyU — (pv) u , Eq. (13. ip and Eq. (|3.2p can be simplified to 



ydxP « dylyr] = dy[y{fidyU - {pv)'u')] (3.3) 



w dy[y{ufldyu - {pv)'h[ + kdyT)]. (3.4) 

To search for a general analogy solution between the momentum and energy equations, 
an extended total enthalpy/temperature, denoted as h^jT^, is defined here to link velocity 
and temperature fields 

he^CpT, = CpT + f{y) — , (3.5) 

where f{y) is called extended recovery function and to be determined later. Using this 
definition and applying the approximations in Eq. (|2.4p . Eq. p.4p can be rearranged to 
(see Appendix B) 



dy[y{{r-l)uf+{mV2)dyf+{l-Vr)qy)] « dymdyh, - {pv)'K)]. (3.6) 



, or integrated from center to to give y{uT — Cp{pv)'T' — qy) ~ 0, where comes from the 
symmetry in the pipe center. This means 



uf ^ Cp{pvyT' + qy. (3.7) 

Similar to the derivation of SRA, a full analogy between the momentum equation 
Eq. (|3.3p and the energy equation Eq. p.6|) can be guaranteed under the constraint of 

dy[y{{f ~ l)uf + {{iu^/2)dyf + (1 - ¥v)qy)] = U^ydxP « U^dy[yf], (3.8) 

where 11^ is a proportional constant with dimension of velocity, and is only associated 
with wall conditions as demonstrated later. To avoid a rigorous proportional assumption 
on fluctuations as that in deriving SRA, we express in terms of u hy — U^u + e , 
where e characterizes the instantaneous deviation of real turbulence from the propor- 
tional relation h^, =^ UwU . Then the energy equation Eq. (j3.6l) becomes 



Un,dy[yT] « dy[y{pdyhe - [pv]' h'^ - {pv)'e')]. (3.9) 

Now we introduced a critical, also unique, assumption: comparing with h^, the devia- 
tion £ is a high-order modification, and hence —{pv)'e' is negligible. Under this assump- 
tion, the momentum equation Eq. p.3p and the energy equation Eq. p.9p can be written 
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in a similar form as 

dy[y{UwT)] dy[y{p,dy{Uwu) -Uw{pv)'u')] (3.10) 
dy[y{Uwf)] dy[y{p,dyhe-Uw{pv)'u')] (3.11) 

Subtracting Eq. (|3.10p from Eq. (|3.11l) and integrating from wall to y, one has y[p.dy{he — 
Uwt)] = const, where const equals zero because of the symmetry at the pipe center. A 
further integration along y gives 

he — UyjU = hyj. (3.12) 

It is worth pointing out the Eq. p.l2l) is a deterministic solution because of the sym- 
metry in pipe, in contrast to the possible solution discussed in SRA for CTBL where no 
symmetry exists. Combining the assumption h^ = U^u + e and Eq. (j3.12[) . we obtain a 
full analogy solution between the momentum equation and the energy equation: 

he — UwU — Uw£ = huj. (3.13) 

Averaging Eq. p.l3p and then subtracting Eq. (|3.12l) . we have e' = 0, which explains 
the reason for using ' as the superscript in e . 

So far, the extended recovery function f(y) is undetermined. Expanding Eq. (|3.12p 
with Eq. p.5|) and calculating its derivative with respect to y, one has: 



CpdyT + fudyu + (u^ /2)dyf=U^dyU. (3.14) 
Applying the equation at the wall yields 



dT 



= -Pr^. (3.15) 



Substituting Eq. (IXTi|) . Eq. ((XT)) and Eq. to Eq. fy can be solved as (see 

Appendix C) 



.(,,.^^-3^, (3.16) 

u u{pv)'u' 

and a full analogy solution is established. 

3.2. Compressible channel flow 
For compressible channel flow, Eq. (|2.2I) and Eq. (12. 3p can be written as 



dxP = dy[Txy - {pv)' u] (3.17) 



= dy[uTxy — {pv)' h[ + vTyy + kdyT]. (3.18) 

Notice that these two equations can be viewed as a consequence of setting y = 1 in 
Eq. p.ip and Eq. (|3.2p : hence all the derivations and conclusions for pipe flow can be 
applied to channel flow by simply replacing y with 1. 

3.3. Compressible turbulent boundary layers at zero-pressure gradient 
For a ZPG CTBL, Eq. and Eq. (gj]) can be written as 



pudxU + pvdyU = dy\{Txy - {pv)' u )\ + dx\Txx - {pu)'u'] (3.19) 
pudxht + pvdyht = dx{uT^) + dx{vT^) + 9^ (mt^) + dy{vT^) 



-dx{pu)'h't - dy{pv)'h't + dx{kdxT) + dy{kdyT). (3.20) 
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Besides those approximates involved in pipe/channel flow, we also neglect the small terms 
associated with streamwise gradient: dx{T^), dx{vTxy), dx{pu)' u , dx{uTxx), dx{pu)' h[, 
dxikdxT). Then Eq. (|3J9l) and Eq. (|3:20l) can be simplified to 



fmdxU + 'pvdyU « dy[p,dyU — {pv)'u'] — dyf (3-21) 

Jmdxht + Jnldyht ~ dy[ufldyU — (pv)' h^] + dy{kdyT). (3.22) 

In CTBL, f can also be treated approximately as a function of only y, and the defini- 
tions for he and Tg, i.e. Eq. p.Sp . still work. Using the approximations in Eq. (|2.4p . the 
energy equation Eq. (j3.22l) can be rewritten as 



pudxhe +pvdyhe W dy{up.dyhe - {pvYh'J + fres (3.23) 

jwhere (see Appendix D) 

fres = pvu^dyf/2 + {f-l)udyf + dy[{PT-l)qy-u^fLdyf/2 + {l-f)uf]. (3.24) 

Similar to pipe/channel fiows, a fully analogy between the momentum equation Eq. (j3.21l) 
and the energy equation Eq. (|3.23|) of ZPG CTBL requires fres = 0, which can be satis- 
fied by solving fijj). However, the exact solution of f{y) is too complex to be determined 
analytically. 

Here we look for an approximate solution of f{y). Using Eq. (I3.14p . fres — can be 
rearranged and decomposed into fres = fresi + fres2 = (see Appendix E), where 



fresi = dy[Cp{pv)'T' - {Uy, - ru){pvyu'] (3.25) 

fres2 = dy[uf - Cp{pv)'T' -%] + [{f - l)u - Uw]dyf 

+{'pv/p)[{Uw - fu)pdyu -f Vvqy]. (3.26) 

Now we introduce a second assumption: comparing with the other terms in Eq. US. 23\) . 
fres2 is negligible. The validity of this assumption is discussed in section |4l Under this 
assumption, the constraint fres = degenerates to fresi — Owhich can be integrated 
from wall to y to give 



Cp{pv) r - {U^ - ru) [pv) u = 0. (3.27) 

This relation is equivalent to the extended recovery function of Eq. p.l6p . So, naturally, 
all the conclusions for pipe and channel flows also work in ZPG CTBL. 

3.4. Adiabatic CTBL with pressure gradient 

For CTBL with pressure gradient, the energy equation takes the same form as that of 
ZPG CTBL. For adiabatic wall^Eq. jXTS]) gives Uyj = 0. Provided that fres2 and dxf-^ 
are still n egligible in thi s flow, he = h^ satisfies Eq. p.23|) independent on the pressure 
gradient (|Gaviglidll987l ). Therefore, the approximate analogy solution for ZPG fiow is 
also suitable for adiabatic CTBL with pressure gradient. 

3.5. General Reynolds analogy 

So far, we obtain full analogy solutions for compressible pipe and channel flows, and ap- 
proximate analogy solutions for ZPG CTBL under arbitrary wall-temperature conditions 
and for adiabatic CTBL with pressure gradient. These solutions have a common form of 
Eq. (|3.13p . To be simple and invariant, we define a general total enthalpy (temperature), 
denoted by hg{Tg), as 



hy = CpTy = CpT + f{y)\ - U^u; f{y) = ^ - (3.28) 

^ u u[pv) u 



12 You-sheng Zhang, Wei-tao Bi, Fazle Hussain, Xin-liang Li and Zhan-su She 



3.0 
2.5 
2.0 
1.5 
1.0' 
0.5 
0.0 
-0.5 



■ CTBLM2.250I.OO 

• CTBLM4.5O01.OO 
A CTBLM4.5O0O.41 
T CTBLM4. 5000.00 

♦ CTBLMe.0001.00 
□ CTBLM1 .500-0.54 
o CTBLM3. 000-0.78 
A CTBLM6.OO0-1.OO 




2.0 
1.5 
1.0 

0.5 
0.0 
-0.5 



-1.0 I ■^^^J ' 1 ' 1 ' 1 ' 1— 1 ' 1 ' 1 ' 1 ' 1 -1.0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Figure 1 . (color online) . Variation of the modified turbulent Prandtl number for compressible 
turbulent boundary layer and compressible channel flow. The DNS of ZPG CTBLs and CCFs 
are simulated by us and Y.Z.W ang (priva te communication) , res pectively. For numerical detail, 
see reference Ref. ljLi et ffl/.ll200lh and Ref. (|Xin-Liang et aZ.|[2006h . 



,where Uw takes the value of Eq. p.l5| ).Then the analogy solutions for compressible wall- 
bounded turbulence can be re-expressed with Tg as 



T — T 



(3.29) 
(3.30) 
(3.31) 

(3.32) 



, where Eq. p.32p is the statistical result of Eq. (|3.31l) , and is because e is a high-order 
perturbation on and u . In this paper, Eq. p.29|) — Eq. p.32p are collectively called 
general Reynolds analogies. 



4. Mean relationships based on the GRA 

In this section, we discuss the relationship between mean velocity and mean tempera- 
ture based on GRA. For convenience's sake, wc first define a modified turbulent Prandtl 
number as 



Prt 



{pv)'u'{dT/dy) 
{pvyT'{du/dy) 



Using this definition and the approximations in Eq. 
Eq. (I3.30p can be expanded to give 



T- 



1 dT 
Prt du 



+ 



dT 
du 



Prt. (4.1) 
(|2.4p . the mean relationship 

f^, (4.2) 



which can be integrated in case of constant . For compressible wall-bounded turbulenc e 
jHuang et all Il995l: [Chosh et all [2OI0I iDuan et all I2OIOI I2OIII: iDuan fc Martini l201l[ ). 
Prt is around one in the wall-normal direction, as shown in Fig. [1] where Prt is plotted 
versus u/ug in contrast to the traditional plot of Prt versus y/S. In the traditional plot, a 
wide plateau of Prt « Prt ~ 0.8 is observed in the outer region of a boundary layer, but 
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corresponding only to a small change of mean velocity. We argue that it is better to plot 
Prt versus u/us because u/us, instead of y /S, is the variable in the temperature-velocity 
relationship, see Ea. (jl.l[) . Only by this plot can the dependence of Prt on the velocity 
variation be better revealed. As shown in Fig. [U Prt is around one up to u/ug w 0.85. 
Hence, an approximation of Prt = const = 1 is used to integrate Eq. (j4.2l) from wall to 
to give 



T 



T 



us_ dT 



fs 



us_ df 



(4.3) 



The derivative of Eq. ()4.3|) with respect to u gives 



dT _ dT 
du du 



us du 

Applying Eq. (14. 4p to Eq. p.28p with the approximation Prt = 1, one obtains 



(4.4) 



r{y) = — , , = -^ij^lw - =j^) 

u u{pv) u u du Prt du 
^Cpdf df f^^fs 2CpdT 
^ u^du^"" du'^^uyi2Cp)^ us du^""'- ^ ' 

Therefore, under the approximation Prt — const = 1, the extended recovery function 
f{y) becomes a constant and this constant only associated with the wall temperature 
condition. Furthermore, Eq. (14. 5p reveals that f{y) consists of two parts: -TT/o!^'' > and 

fs" • Similar to the recovery factor r defined for an adiabatic wall, the first part 
can be understood as a nominal recovery factor at a non-adiabatic isothermal wall. The 
second part denotes the contribution to f{y) from the non-zero heat flux at the wall, 
which becomes important if the wall is over-heated or over-cooled, such as that in CCF or 
CPF. For an adiabatic wall, f naturally reduces to the classical recovery coefficient, saying 
f = r = {fr — fs)/ {us /2Cp). Owing to these reasons, f is called extended recovery factor, 
instead of function, hereinafter. Correspondingly, an extended recovery temperature can 
be defined as ff = fs + rus/ {2Cp). With this definition, Ea. (j4.3p can be rewritten as 
Equation Section (Next) 

rrri rri rji — rji rri / — \ ^ 

-L ijf} -L -L uj Ui J- (5 J. I Ui 



^ ^ . ^ - ^ ^ ,^ ,T,^T + f—, (4.6) 
Is Is J^s us Is \us J 2Cp 

which has the same form as Crocco-Busemann relation and Walz's equation except the 
different value of f. For an adiabatic wall, Eq. (14. 6p reduces to the widely accepted Walz's 
equation. This general relationship shows that the mean temperature has a quadratic 
dependence on the mean streamwise velocity at arbitrary wall-temperature conditions, 
and the effect of wall temperature is only to change the value of f. Note that Eq. (j4.6p is 
universal for CCF, CPF and CTBL. Figl2]and Fig E] compare Eq. (jM]) with DNS data 
of CCF, CPF and CTBL, which show a collapse between the current theory and the DNS 
data. In Fig. [51 Eq. (|4.6p displays a much better performance than Walz's equation for 
non-adiabatic CTBL. 

f in Eq. ()4.6p requires exact wall information to be calculated, which is inconvenient in 
the engineering application. Therefore, we investigate the influence of wall-temperature 
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Figure 2. (color online). Mach number (left) and wall-temperature effects (right) of the 
temperature-velocity relationship Eq. (|4.6|l in CTBL. 
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Figure 3. (color online). Comparison between the temperature- velocity relationship Eq. (|4.6[l 
and DNS data in compressible channel flows (left) and compressible pipe flows at different 
adiabatic d egrees and Mach numbers. The DNS of CCFMl.^QQ — 0.71 com es from reference 
l|Huang et a/.iil995.) . and CPFMl.ZQQ - 0.76 comes from reference l|Ghosh et ffl/.„20ia) . 



on r. To do this, we introduce the vifell-known Reynolds analogy factor defined as 

Qyw^S 



(4.7) 



where C 



f 



is skin- friction coefficient, Ch 



is heat-transfer coef- 



ficient, i.e. Stanton number (|Smits fc Dussaugell2006l ). Using this definition, the general 



recovery factor r can be rewritten as (see Appendix F) 

f = r[sPr+(l - sPr)e]. 



(4.8) 



An outstanding property of Eq. (|4.8p is that all the param eters in the right hand side of 



Eg. (14.81), except 6, are alm ost constant. In the literature (|Bradshawlll977t iDuan et 
20ld : Duan fc Martm 2011 ). the Reynolds analogy factor s varies between 1.0 and 1.2 



without a clear trend with respect to T^,, i?e, M, chemical reaction and enthalpy condi- 
tion. Here, we use sPr = 81 ^ J^, which is set by averaging the DNS data of us and 
several others ( Duan et al. 2010l : iDuan fc Martini 1201 ll : fchosh et a/.l l2010': Hu ang et al\ 
as shown in Fig.(|4]). Using Eq. ()4.8|) and the definition of 6, Eq. (|4.6p can be 
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rearranged to (see Appendix G) 



T 



fs 



+ 



T —T 



us Ts 



us 



fi- 



us 



(1 



us 



«(- 



us 



(4.9) 



which is consistent with Duan et.al's empirical relation Eq. (|2.29p except the slight dif- 
ference between a = sPr ^ 0.81 in EqMM and a = 0.8259 in Eg. I^Mi){,see fig.Q). 



Using Duan et.al (jPuan fc Martini l201l[ )'s (jPuan fc Mart i'nll20 111 ) definition about T*^^, 
Eq. ()2.29|) can further be rearranged to: 



y us 



u 

US' 



u 

US' 



s Pr 



0.81. 



(4.10) 



In this form, the T*y is a deterministic function of u/us- Hence, a collapse of different 
DNS data is expected. This is confirmed by figlH where the non-adiabatic DNS data of 
CPF, CCF and CTBL are plotted together. Indeed, an excellent collapse of DNS data 
to Ea. (HJU|) is observed in the range of M = 1.26 ~ 9.4 and 9 = 0.6 ~ -1.0. Since 
the DNS by Duan et.al also considered the effect of real gas (enthalpy) and surface 
catalysis, such a collapse shows that the validity of Eq. (|4.6p is independent on T^, 
Re, M, flow systems s urface catalysis and enthalpy condition. Noting that Duan et al. 
( Duan fc Martin 1201 ll) already observed that the temperature- velocity relationship was 
influenced by both enthalpy condition and surface catalysis, here we clarify that this 
influence acts by changing the wall temperature (thus f ), without involving any more 
complex reason. Additionally, we point out that the collapse in Fig|3|can only be realized 
by setting r ss 0.9 for CPF and CCF, which means that wall-bounded turbulence shares 
the same adiabatic recovery factor. 



5. Fluctuation relationships based on the GRA 

From Eq.(|32i) and Eq.dSjnOj), we have 





u{pv)'u' 



T 



(5.1) 
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-^U-? ^-f. (5.2) 
u{pv)'u' J 2 du ^ ' 

Subtracting Eq. (j5.2l) from Ea. (|5.1l) and using the approximations in Eq. (j2.4l) . we obtain 




T 



(pvVT' 
{pv)'u' 



If the contribution from e to the r.m.s. of T is neghgible, Eq. (|5.3I) yields 



T'2 _^ ^{pv)'T' _ ^pv'T 

u 2 {pvyw 



pv'u' 



(5.3) 



(5.4) 



where subtraction is for CTBL and plus for CPF/CCF, a convention in this subsec- 
tion. The different signs are owing to the different temperature distribution along the 
wall- normal direction between CTBL and CPF/CCF. In CTBL, the high-temperature 
region is located near the wall, where the correlations pv'T' and pv'u' have opposite 
signs because of the ejections of lo w-speed/high-ten i perature fluids and the swee ps of 
high-speed/low-temperature fluids ( Duan et aLl 12010 : IPirozzoli fc Bernardinil 120111 ). An 
opposite situation happens in CCF and CPF wh ere the low-temperat ure region is located 
near the wall. As pointed out by Guarini et.al ( Guarini et al. 200Ci ). Eq. (|5.4l) expresses 
an analogy between the rates of turbulent heat and momentum transfer. If we further 
introduce the approximations of {pv)'T' « pv' T' and (pv)'u ' » pv 'u , Eq. (|5.4p yields 
Rv'u' ~ iRv'T', which is the basis of HSRA ( Guarini et al. 2000). The equal correla- 
tions reveal that the turbulence- induced wall- normal transfer of heat is closely associated 
with that of momentum. In other wor ds, the tempera ture fleld is a passive scalar field, 
consistent with Morkovin's viewpoint (|Morkovinlll962r ). 

Combining Eq. (|2.4p and Eq. (|5.4p . a relationship between and can be obtained: 



1 



(5.5) 



{j~l)Ma^V7^/u ±Prt{l-dTt/dT) 

This relation has the same form as that of HSRA, except replacing c = Prt in Eq. (|1.2p 
with c = Prt- Here, we call Ea. (l5.5p the modified HSRA (MHSRA). A comparison 
between DNS and Eq. (|5.5l) is demonstrated for CTBL in FiglS] To visualize clearly, FiglS] 
plots the ratio of the left-hand side of SRA, HSRA and MHSRA to their right-hand side 
at different M and 8. Fig (S] shows that MHSRA is in excellent agreement with the DNS 
data throughout the boundary layer, while SRA is valid only to adiabatic wall and below 
0.76. For SRA, an increased deviation from unity is observed in flows with colder walls. 
HSRA is plotted in the inset of the right figure. It is shown that MHSRA has a slight 
improvement over HSRA. 

Using Eq.f 15. 5p . we can explain the validity of Eg. (12. 181) at adiabatic wall, which is 
said to satisfy under an incorrect assumption as mentioned in sec l2.2l With Eq.( 14. 5p . a 
simple calculation gives 



1 



±PTt{dTt/dT -1) 



Cp 1 dT 
Prt u du 



CpdT 
u du 



(5.6) 



For quasi-adiabatic wall, the right hand side of this relation is r « 0.9, close to the unity 
value of SRA, and, subsequently, leading to the good performance of SRA. For non- 

adiabatic wall, however, the deviation from unity, caused by local modification 
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Figure 5. (color online). A comparison between SRA (left) and MHSRA (right) at different 
M and 6. The inset in the right figure is HSRA. 



and global change of f , becomes more and more significant with the increase of heat flux 
at wall, resulting in the failure of SRA. 

For the correlation coefficient |i?„'T'|, on one hand, the existence of high-order per- 
turbation e in the instantaneous fluctuation relation Eq. (|5.3|) excludes the possibility 
of |i?„'T'| = 1, an incorrect conclusion of SRA,ESRA,GSRA and HSRA; on the other 
hand, it prevents a simple formulation about To avoid this dilemma, we recur 

to MHSRA. Simplifying the left-hand-side of Eq.(l53|), we have (?2/Cp)(v^/v^) = 
±Pr7((9rt/ar - 1). Substitution of this relation into Eq.^^ yields 



±- 



1 



t;^/t"^ 



1 



PrtidTt/dT-l) 



Prt{dft/df -I) 



(5.7) 



which gives a connection between Ru'T' and Prt- Theoretically speaking, Ea. (|5.7|) has 
the same accuracy as that of MHSRA. 



6. Discussion 

Now we discuss the key approximations and assumptions involved in GRA. For CCF 
and CPF, GRA is a full analogy solution if the unique assumption of negligible e is valid. 
For CTBL, additional assumptions of negligible fres2 = and streamwise derivatives are 
applied to give an approximate analogy solution. The assumption about e is indirectly 
confirmed by the excellent agreement between DNS and GRA. The approximation of 
negligible streamwise derivatives is widely adopted in the study of turbulent boundary 
layers, which requires that the flows are slowly developing in the streamwise direction. 
The requirement prohibits an incautious application of GRA in CTBLs that have strong 
acceleration or temperature gradient along the wall, large local wall curvature, intense 
interactions with an impinging shock or a protruding object, etc. The self-consistence 
of approximation fres2 ~ can be verified. Under the solution of GRA, the residual 
error fres2 is equivalent to fres2 — T'dy[{^ ~ + {Pt: —\)dy{qy) (see Appendix H). 
Approximation fres2 ~ is equivalent to require that fres2 is a small term compared 
with the dominant terms of Eq. (|3.23|) . which is supported by DNS data (see Fig[S]). Note 
fres2 is negligible above the near- wall region, where fdyU and dy{qy) is not important 
and f ^ Pr. A possible deviation of GRA from DNS may observed in the near wall region 
where a small residual of fres2 exist (see Figl6|). This deviation increases for cold walls 
that have large dy{qy) near the wall, as shown in Fig.O However, the difference between 



18 You-sheng Zhang, Wei-tao Bi, Fazle Hussain, Xin-liang Li and Zhan-su She 

0.3 
0.2 
0.1 



0.0 
-0.1 
-0.2 

-0.3 

1E-3 0.01 y/g 0.1 1 

Figure 6. (color online). Distributions of the budgets of Ea. (|3.23p along the wall- normal 
direction. The budegts are renormalized by the free-stream value. 

MHSRA and DNS in the near wall region is not observed to the mean temperature- 
velocity relationship of GRA. Finally, for the mean relationship, the current theory is 
independent on Pr but requires Prt — const, and the approximation of unity Prt gives 
a rather accurate mean relationship as shown in SecUl It is worthy investigating whether 
or not the assumption Prt — const works in other flow systems. 

GRA demonstrates universality to various compressible wall-bounded turbulent flows 
under various wall temperature conditions. For the mean temperature-velocity relation- 
ship, GRA shows that the coupling between the mean velocity and mean thermal fields 
is same for various flows: the general total enthalpy (temperature) hg{Tg) keeps the con- 
stant value of hgwiT^) along the wall-normal direction. hg{Tg) successfully integrates 
the effects of viscous dissipation and heat flux through the wall on the energy transfer 
between kinetic energy and fluid enthalpy in wall-bounded flows. For the fluctuation 
relationship, the unique difference between internal flow (CCF/CPF) and external flow 
(CTBL) is the sign in Ea. (|5.5p . As discussed in seclSl the difference is caused essentially 
by the different distributions of temperature, a passive scalar in nature determined by 
boundary conditions. 

The universality of GRA needs a further discussion. We argue that the reason behind is 
the common turbulence mechanism in wall-bounded turbulence. Specifically, the essential 
dynamics of large scale eddy movements or coherent structures in wall-bounded turbulent 
flows are similar, without significantly infiuenced by compressibility, wall-temperature 
condition and flow situation far away from the wall. The turbulence transport of mo- 
mentum and heat is dominated by the invariant vortical structures; thus a universal 
analogy solution for the velocity and temperature fields, formally like Walz's equation 
that is valid only to adiabatic CTBL and Duan's relation that applies to many more 
flows, is promised. GRA seems to be such a universal solution. 

The adiabatic degree Q introduced in GRA is a well-founded dimensionless param- 
eter for characterizing the wall-temperature effects in non-adiabatic compressible wall- 
bounded turbulent flows. In these flows, three control parameters are generally involved: 
M, i?e, and a dimensionless wall temperature that, we argue, has not been well defined 
up to now. To study the M-effect, which is an overwhelming topic in the literature of 
compressible turbulence, one has to compare different M flows at a same Re and a same 
dimensionless wall temperature, whose definitions are not trivial owing to the nonlinear 
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couplings between the different physical mechanisms they represent. We recently pro- 
posed a definition of Res^^ and showed it was appropriate for studying the M effects in 
adiabatic CTBLs. Note that the adiabatic degree 9 in Eq. (|4.8p is the unique parameter 
determining f, and hence determining the temperature- velocity relationship of Eq. (I4.6p . 
So Q poses as a well-founded dimensionless parameter that can be applied to study the 
M-efFects of non-adiabatic compressible wall-bounded turbulence. 

The current theory can be applied to make quantitative predictions of the mean fields 
of compressible wall-bounded turbulence. A long lasting wish is to quantify the com- 
pressible wall-bounded turbulence with the much more fruitful results of incompressible 
wall-bounded turbulence, such as the recent formulations obtained by She et al. on i n- 
comp ressible channel flow, pipe flow and turbulent boundary layer ( She et al. 201lQl f6. 
20101) . The key question is to quantify the Al and 9 effects. For an adiabatic compressible 



wall-bounded turbulent fiow (i.e. 9 = 1), the influence of M can be accounted by con- 
sidering the variation of the mean fluid properties, which is the well-known Morkovin's 
hypothesis. Recently, we specifled this hypothesis by considering both the mean density 
and viscosity and successfully removed the M-effects of a series of important flow quanti- 
ties (.You-sheng et a/.,. 2012.) . Because of the nonlinear coupling between the velocity and 
thermal fields, these fiow quantities cannot be predicted by those of the corresponding 
incompressible flows without the information of the mean fluid properties. GRA docu- 
mented in this paper provides a solution for this problem. Now, combining GRA and the 
Ajf-invariants rescaled by the mean thermal quantities, we are able to predict the mean 
profiles of velocity, temperature, density, as well as the Reynolds stress of compressible 
wall-bounded turbulence by using the formulations in the corresponding incompressible 
fiows. We will report the results elsewhere. 



7. Conclusions 

In summary, we presented a general Reynolds analogy (GRA) between the Reynolds 
averaged momentum equation and energy equation of compressible wall-bounded turbu- 
lent fiows. A full analogy is obtained for CPF and CCF, and an approximate analogy is 
obtained for CTBL. Several relations between the temperature and velocity fields are de- 
rived by GRA and validated by DNS data of CPF, CCF and CTBL. It is shown that the 
GRA is superior to the previous analogy theories for its higher accuracy, invariance with 
wall temperature condition and fiow system, and solid physical basis. The GRA opens a 
door to the long lasting wish of predicting the mean fields of compressible wall-bounded 
turbulence with the information of the corresponding incompressible fiow. The method 
to derive GRA may also be suitable for the study of other fiow systems associated with 
thermal fields, such Rayleigh-Benard fiow. 

We thank J. Chen, X. Chen, Y. Wu, Y.Z. Wang, H. Liu, D.X. Fu, Y.W. Ma, X. Liang 
and Z.W. He for helpful discussions, SSC(Shanghai) and NSCC(Tianjin) for providing 
computational resource. This work is supported by NSFC under Grant 90716008 and 
10921202, by MOST 973 Project 2009CB72410, and by CAS Program KJCX2-EW-J01. 



Appendix A. 

Reference ( Ghosh et al. 2010l ) gives Navier-Stokes equations in cylindrical coordinates, 
of which continue equation, streamwise momentum equation and energy equation can be 
averaged in time, axial and azimuthal direction to give, respectively. 



dyiypv] = 



(Al) 
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dy [ij'pmj] = -ydxP + dy [j/t^] (A 2) 

dy [yphtv] = dy [ykdyT + Wyy + uTxy] (A 3) 

Continue equation means JrU = 0, and then one has pvf = JrUf + (pv)' f = (pv)' f . 
Substitution of this relation into momentum and energy equation gives immediately 
equation Eg. dXTj) and Eq. ((3?2|) . 

Appendix B. 

Using the definition about h^ and approximation Eq. (j2.4l) . one has 

ujldyU w fldy[he — (f — l)i2^/2 — CpT] 

= jldyihe) - (f - l)uifidyu) - p,{u^/2)dyf - Prkdyf (B 1) 

{pv)'ht = ipv)'{he - (f - 1)mV2)' w {pv)'h'^ - (f - l)u{pvyu (B 2) 

Substituting the two relation into Eq. p.2p gives 

O^dy [y{fldy{he) - (f - l)uf - fl{u^ /2)dyf - Pr kdyT - {pv)' h'^ + kd^)] (B 3) 
Rearrangement last relation gives Eq. p.Gp 

Appendix C. 

Integration the constraint Eq. (|3.8p from center to y gives 

y[{f - l)uf + {pu'^/2)dyf + (1 - Vi)qy - U^f] w co7ist (C 1) 

In the center, the symmetry imposes qys ~ 0, {dyf)s — and ~ 0, and naturally 
const — 0. Then Eq. (jC 1|) becomes 

(f - l)uf + {iiu^/2)dyr + (1 - VT)qy -U^f^Q (C 2) 

After multiplying Eq. (|3.14|) by p., one has 

{pv? /2)dyr — {Uw — fn)pdyV + Prg^ (C 3) 

Substituting this relation and Eq. (|3.7p into Eq. p.8|) yields 

(f - l)[Cp{pvyT' + qy] + ({/„, - fu)pdyv + PTqy + {l- Pr)^j, - C/„f w (C4) 

Rearrangement of last equation with the definition about f gives 

(f - l)[Cp{pvyT'] + f{qy - vf) + {U^ - fu){pvyu' « (C 5) 

With Eq. p.7|) . one has qy — uf — —Cp{pvyT' , substitution of which into last equation 
gives 

Cp{pvyT' ~ {Uuj - ru){pvyu K. 0. (C 6) 

Solving this equation in terms of f gives Eq. (|3.16p 

Appendix D. 

Using the definition about he, approximation Eq. ip^ . Eq. (|B II) and Ea. (jB2|) . Ea. p.22p 
becomes 

-pndx [K - (f - l)uV2] + pvdy [K - (r - 1)^2/2] w 

dy[fldy{he)-{f-l)uf-p{u'^/2)dyf-FTkdyT-{pvyh'e]+dy{kdyT) ■ ^ 
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Rearrangement of last equation gives 



pudxhe +pvdyhe « dy[fidy{he) - {pvyh'^] + fres, (D 2) 

where 

fres = (f — l)upudxU + (f — l)'pvudyU + 'pvu^dyf /2 

+dy[{PT -l)qy - u^fldyf/2 + f)uf]. (D3) 

Note that the streamwise derivative 'pudxU can be written with normal derivative by 
momentum equation, i.e. JmdxU ~ dyf — 'pvdyU. Then can be further written as 

fres = [{r - l)u\[dyf - pvdyu] + [{f - l)u]pvdyu + [pvu'^/2]dyf 

+ dy[iPT-l)qy - fl{uy2)dyf + {1 - f)uf]. (D4) 

The cancellation of [(f — l)u]'pvdyU in last equation gives Eq. p.24[) . 

Appendix E. 

Rearrangement of Eq. p.24p gives 

fres = {pv/fl)iu^fldyf/2) + dy[iPT-l)qy - {u^fidyf/2)] + fdyHl " f)u]. (E l) 

After replacing u^fldyf/2 with Ea. (|C 3p . a further rearrangement yields 



fres = {pv/lJ')[{Uw - ru)fidyU + Pr Qy] + dy[uT - Qy - {Uw - ru){pv)'u'] 

+ [{f-l)u-U^]dyf. (E2) 

Therefore, one has 



fres = dy[Cp{pv)'T' - {Uui - ru){pv)'u'] + [{r ~ l)u - Uj\dy 



T 



+dy[uf - CpipvYT' -q^] + {pv/fL)[{U^ - fu)pdyu + Pi qy\. (E 3) 
A decomposition gives fresl and fres2 



Appendix F. 

With Ea. (|4.5p . the f can be further written as 

Tr -fs ff^-fs , Us df 



Because of jr^fr |? 



u]/{2Cp) \ fr-fs ' Tr-Ts du 
Tw —Ts — us dT 

Tj. — Ts Tj. — Ts Tj- — Tyj du 

= sPr, = 1 - and = 6, one has 



(Fl) 



f{y) = r[e + (1 - e)sPr] = r[s Pr+(1 - sPr)e] (F2) 

Appendix G. 

With the definition about Tf and Tr, one has 

Tr=Tr + {f/r - \){rul/{2Cp)) ^ Z + (.sPr-l)(l - e)(T, - fs) 

= T, + (sPr-l)(T, -T„) (Gl) 
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Substituting this relation into Eq. ()4.6|) , one has 



T ^ Tr -T^: + {sFr-l)(Tr ~T^) u ^ Tg - Tr - [sVy -l){Tr - T^) ^^^^ 



Ts Ts Ts us Ts us 

Rearrangement of last equation gives 

2 



T T T — T 



Ts Ts Ts 
jwhich gives Eq. 



sPr^ + (l-sPr) ( — 
us \us 



Ts — Tr ( u ^ ^ 



Ts \us 



(G3) 



Appendix H. 

Using the relation Eq. (|3.27p . Eg. (13.261) can be rewritten as 

/res2 = dy\uf ~ {Un, - fu){fldyU - t) - qy] + [(f - l)u - Uj\dyf 

- fu)fldyU + Prqy] (H 1) 

Rearrangement of this equation gives 

fres2 = Tdy[{l - f)u] + (Pr -l)dy{qy) + dy[{Un, - fu)jidyU + Prgj^] 

+ipv/fl) [{Uuj - fu)fldyU + Pr qy] (H 2) 

With the fourth term of Eq. ()4.5p . one has {Uw — fu)pdyU + Pr qy = 0. Then /, es2 can be 
further simplified to 

fres2 = fdy[{l - f)u] + {Pv -l)dy{qy) (H3) 
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